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Use  of  a  Two-Dimensional  Flow  Model  to  Quantify  Aquatic  Habitat 

by 

D.  Michael  Gee1  and  Daniel  B.  Wilcox2 


Introduction 

^This  paper  describes  the  use  of  a  numerical  two-dimensional  flow 
model  to  evaluate- the  impacts  of  potential  hydropower  retrofits  on 
downstream  flow  distributions  at  Lock  and  Dam  No.  8  XFig.  1)’ on  the 
upper  Mississippi  River.  The  model  used  (RMA-2  <6H  solves  the 
complete  Reynolds  equations  for  two-dimensional  free-surface  flow  in 
the  horizontal  plane  using  a  finite  element  solution  scheme.  RMA-2 
has  been  in  continuing  use  and  development  at  the  Hydrologic 
Engineering  Center  and  elsewhere  for  the  past  decade .  WyBrd.S)* 
Although  designed  primarily  for  the  simulation  of  hydraulic 
conditions,  RMA-2  may  be  used  in  conjunction  with  related  numerical 
models  to  simulate  sediment  transport  and  water  quality  (5,7)/  In 
this  study,  velocity  distributions  were  evaluated  with  regard  to 
environmental,  navigational  and  small-boat  safety  considerations. 
Aquatic  habitat  was  defined  by  depth,  substrate  type  and  current 
velocity.  Habitat  types  were  quantified  by  measuring  the  areas 
between  calculated  contours  of  velocity  magnitude  (isotachs)  for 
existing  and  project  conditions.  The  capability  for  computing  and 
displaying  isotachs  for  the  depth-average  velocity,  velocity  one 
foot  from  the  bottom  and  near  the  water  surface  was  developed  for 
this  study.  The  product  of  this  study  effort  is  an  application  of 
the  RMA-2  model  that  allows  prediction  of  structural  aquatic  habitat 
in  hydraulically  complex  locations.  Elements  of  the  instream  flow 
group  methodology  (1)  could  be  incorporated  to  provide  detailed 
predictions  of  Impacts  to  habitat  quality. 

Calibration  of  the  numerical  model  to  field  measurements  of 
velocity  magnitude  and  direction  is  also  described.^ 

Project  Description; 

Lock  and  Dam  No.  8  is  located  on  the  upper  Mississippi  River 
679.2  river  miles  (1094  km)  above  the  mouth  of  the  Ohio  River,  23.3 
miles  (37.5km)  below  Lock  and  Dam  No.  7,  and  31.3  river  miles 
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(50. Aka)  above  Lock  and  Dan  No.  9.  The  lock  li  on  the  left  bank,  or 
Wisconsin  side  of  the  river  0.3  miles  (0.5  ka)  below  the  village  of 
Genoa,  Wisconsin,  which  is  on  State  Highway  No.  35,  17  wiles  (27  ka) 
south  of  the  City  of  La  Crosse,  Wisconsin. 

The  aain  lock  is  110  feet  (33.6  a)  wide  and  600  feet  (182  a) 
long,  and  the  upper  gate  bay  of  an  auxiliary  lock  is  provided  in  the 
event  it  becoaes  necessary  to  add  another  lock  in  the  future.  From 
the  river  wall  of  the  auxiliary  lock,  a  movable  daa  section 
consisting  of  five  roller  gates,  20  feet  (6.1  a)  x  80  feet  (24.4  a), 
and  10  tainter  gates,  IS  feet  (4.6  a)  x  35  feet  (10.7  a),  extends 
across  the  swin  channel  to  the  right  bank  of  the  river. 

Hydropower  installations  were  being  considered  by  the  Corps  of 
Engineers  for  the  storage  yard  at  the  west  abutaent  or  in  the 
auxiliary  lock  location.  In  addition  to  the  Corps  plans,  privately 
funded  installation  of  floating  powerplants,  "hydrobarges",  is  being 
considered  for  soae  of  the  roller  gate  bays. 


A  key  coaponent  of  the  successful  application  of  RMA-2  is  the 
developaent  of  a  well  designed  finite  element  network  that  accurately 
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represents  flow  boundaries,  bottom  topography,  inflow  and  outflow 
locations,  and  contains  enough  detail  to  resolve  flow  patterns  of 
interest.  A  study  area  was  selected  extending  from  the  downstream 
face  of  the  dam  downstream  approximately  2.2  miles  (3.S  km).  The 
downstream  boundary  was  set  at  approximately  river  mile  677  (km 
1090).  This  location  was  judged  to  be  sufficiently  far  downstream 
from  the  dam  that  differing  flow  distributions  at  the  structure, 
representing  the  different  project  configurations,  would  not  affect 
the  water  surface  elevation  (assumed  horizontal)  at  the  boundary. 

Bathymetry  for  the  study  reach  was  provided  at  1:2400  map  scale 
with  a  5  ft.  (l.S  m)  contour  interval.  A  finite  element  network  was 
constructed  to  define  the  study  area  with  adequate  detail  at  the  dam 
face  to  distinguish  releases  from  individual  spillway  bays.  Several 
wing  dams  and  training  structures  within  the  study  reach  are 
reflected  in  the  network  as  well  as  a  large  island  in  the  lower  half 
of  the  reach  (Fig.  2).  Curved-sided  elements  were  used  along  lateral 
flow  boundaries  to  allow  for  smooth  tangential  flows  at  the  sides 
(slip  flow  boundary  conditions).  The  network  was  revised  soawwhat 
during  calibration.  The  final  network  consisted  of  37S  elements  and 
1,189  nodes.  Both  triangular  and  quadrilateral  elements  were  used. 

Boundary  Conditions: 

Boundary  condition  selection  was  straightforward  for  this  study. 
Releases  from  the  spillway  bays  and/or  powerplants  were  described  as 
flows.  The  downstream  boundary  condition  was  a  specified  water 
surface  elevation,  assumed  horizontal  across  the  river.  That 
elevation  was  varied,  depending  upon  the  flow  magnitude,  based  upon 
information  provided  by  the  St.  Paul  District.  Lateral  boundaries 
were  generally  "slip”  conditions  which  allow  flow  parallel  to  the 
boundary.  In  some  cases  where  sharp  corners  exist  on  the  boundary  a 
zero  velocity  was  specified  (stagnation  point  boundary  condition). 

All  simulations  assumed  steady  flow,  although  RMA-2  can  readily  be 
used  to  simulate  unsteady  flow  problem  as  well. 

Calibration: 

An  essential  ingredient  to  any  modeling  effort,  be  it  physical 
or  numerical,  is  calibration  of  various  model  parameters  or 
coefficients  such  that  model  results  acceptably  reproduce  prototype 
observations.  The  parameters  available  for  calibration  of  RMA-2  are 
bottom  roughness  and  turbulent  exchange  coefficients.  Calibration 
strategy  typically  consists  of  estimation  of  bottom  roughness  values 
(based  upon  knowledge  of  bed  material)  followed  by  adjustsmnt  to 
reproduce  prototype  head  loss  within  the  study  reach.  Turbulent 
exchange  coefficients  are  then  adjusted  to  reproduce  details  of 
prototype  velocity  distributions.  Thus,  it  is  desirable  to  have 
prototype  measurements  of  velocity  (magnitude  and  direction). 


Velocity  measurements  were  obtained  at  ten  transects  in  the 
study  reach.  Coverage  was  very  good  consisting  of  up  to  five  point 
measurements  at  each  of  up  to  ten  locations  across  each  transect. 
This  calibration  data  set  was  obtained  during  the  period  15-16 
August  1984  at  a  river  flow  of  25,400  cfs  (711  m’/sec)  using  a 


Figure  2 

Finite  Element  Network  -  Flow  is  from  Top  to  Bottom 


Price  Meter.  Although  the  Price  meter  only  yields  directional 
information  for  near-surface  measurements,  these  directions  were 
assumed  to  apply  to  the  entire  water  column  at  each  point. 

Based  upon  these  velocity  smasurements ,  roughness  coefficients 
were  adjusted  to  match  the  flow  distribution  around  the  island 
indicated  by  the  data;  approximately  65%  to  the  west,  35%  to  the 
east.  Additionally,  the  transverse  velocity  distributions  at 
certain  transects  were  not  symmetric,  with  lower  velocities  on  the 
right  (west)  side  than  on  the  left.  This  did  not  appear  to  be 
entirely  due  to  the  presence  of  the  right  bank  wing  dams ,  therefore, 
the  bottom  roughness  was  increased  along  the  right  side  of  the 
channel  from  the  wing  dams  to  about  the  head  of  the  island  to 
reproduce  the  skewed  velocity  distributions  (see  Figure  3). 

Roughness  used,  in  terms  of  Manning  n-values,  varied  from  0.030  to 
0.060. 


The  turbulent  exchange  coefficients  were  set  at  25  pound-second/ 
feet3  (120  kg-sec/m1 2)  for  small  elements  and  50  pound-second/ 
feet2  (240  kg-sec/m2)  for  all  others.  This  set  of  coefficients 
and  network  yielded  well-behaved  numerical  solutions  that  generally 
converged  in  four  iterations.  Internal  continuity  checks,  an 
indicator  of  network  adequacy,  were  generally  within  +  5%. 

Calibration  of  a  numerical  two-dimensional  flow  model  such  as 
RMA-2  is  a  judgmental  process.  That  is,  no  objective,  or  statisti¬ 
cal,  measures  of  error  between  measured  and  simulated  velocities  or 
water  surface  elevations  are  used.  In  this  study,  plots  of  velocity 
vectors  (observed  and  simulated)  were  examined  visually. 

Coefficients  were  then  adjusted  until,  in  the  judgment  of  the 
modeler,  an  acceptable  balance  was  achieved  between  reproduction  of 
the  measurements  and  physically  realistic  coefficient  values.  Some 
network  amdif ications  were  awde.  Additional  detail  was  added  in  some 
areas  to  better  reproduce  observed  velocity  gradients.  A  portion  of 
the  simulated  and  measured  flow  fields  is  shown  in  Figure  3. 

Production  Simulations; 


All  project  (i.e.,  hydropower)  siswlations  were  made  in 
conjunction  with  a  corresponding  existing  (i.e.,  no  hydropower) 
condition  simulation.  Thus,  project  impacts  are  defined  as 
differences  between  project  and  existing  conditions.  Approxiswtely 
25  different  conditions  were  analysed. 

Results  of  the  numerical  simulations  were  presented  in  the 
following  fonnats  for  further  analysis: 

1.  Plots  of  velocity  vectors  approximately  one  foot  (0.3  m) 
above  the  bottom. 

2.  Plots  of  velocity  vectors  near  the  water  surface. 

3.  Plots  of  velocity  vectors  of  vertically  averaged  velocity 
(the  model  output) . 
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4.  Isotachs  (i.e.,  contours  of  equal  velocity  magnitude)  one 
foot  (0.3  m)  above  the  bottom. 

5.  Isotachs  near  the  surface. 

6.  Pathline  plots  depicting  paths  taken  by  particles  traveling  at 
the  vertically-averaged  velocity  released  at  various  points. 

These  plots  are  not  reproduced  herein  in  their  entirety  due  to 
space  limitations;  examples  of  the  above  graphical  displays  are 
shown  in  Figures  4-6. 

The  velocities  near  the  bottom  and  surface  were  estimated  from 
the  computed  vertically-averaged  velocity  by  fitting  a  classic 
logarithmic  velocity  profile  at  each  point. 


Interpretation  of  Results 


Kiverine  habitat  types,  defined  by  substrate  type,  depth,  and 
current  velocity,  vary  in  areal  extent  and  distribution  with  river 
discharge.  The  location  and  extent  of  habitats  in  the  tailwater  of 
lock  and  dam  8  were  quantified  using  the  model  simulations  of  current 
velocity  at  one  foot  above  the  bottom.  Model  results  indicate  that 
hydropower  operation  would  affect  the  pattern  of  habitat  types  only 
within  2,000  feet  (610  m)  of  the  dam.  The  greatest  net  change  in 
aquatic  habitat  between  existing  and  with-project  conditions  would 
occur  at  the  lowest  river  discharges.  Model  results  indicate  that 
hydropower  operation  would  cause  little  disruption  of  habitat  types 
and  fish  distribution  during  the  high  discharge  period  in  spring 
when  saugers,  walleyes,  and  anglers  congregate  below  the  dam. 


An  analysis  similar  to  that  used  to  quantify  aquatic  habitat  at 
lock  and  dam  8  could  be  used  to  evaluate  effects  of  hydropower 
operation  on  navigation,  by  employing  simulations  of  surface 
currents . 


Conclusions 


The  application  of  KMA-2  described  herein  represents  a  new  use 
of  the  model  that  had  previously  been  used  for  hydraulic,  water 
quality,  and  sediment  transport  studies.  The  importance  of  graphics 
post  processors  for  interpretation  of  the  results  of  two-dimensional 
flow  simulations  is  highlighted  in  this  study.  The  methods  of 
extracting  meaningful  information  from  simulated  flow  fields 
illustrated  herein  should  prove  valuable  for  habitat  evaluation  in 
general . 
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Figure  4  Example  Vector  Plot 

Spillway  Discharge  ■  3000  cfs  (84  cms) 
Powerplant  Discharge  -  7000  cfs  (196  cms) 
(located  in  auxiliary  lock) 
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Center  in  Maintaining  Widely  Used 
Hydrologic  and  Water  Resource 
Computer  Models,  B.  S.  Eichert, 

Nov  78,  16  pp. 

TP-57  Flood  Damage  Assessments  Using  Spatial  ADA-106  256 

Data  Management  Techniques,  D.  W.  Davis 
and  R.  P.  Webb,  May  78,  27  pp. 

TP-58  A  Model  for  Evaluating  Runoff-Quality  in  ADA-106  257 

Metropolitan  Master  Planning, 

L.  A.  Roesner,  H.  M.  Nichandros, 

R.  P.  Shubinski,  A.  D.  Feldman, 

J.  W.  Abbott,  and  A.  0.  Friedland, 

Apr  72,  81  pp. 


TECHNICAL  PAPERS  (TP) (Continued) 


HEC 

NUMBER 


TITLE 


HEC 

PRICE 


NTIS 

NUMBER 


E2 .00  Etch 


TP-59  Testing  of  Several  Runoff  Models  on  an 
Urban  Watershed,  J.  Abbott, 

Oct  78,  53  pp. 

TP-60  Operational  Simulation  of  a  Reservoir 
System  with  Pumped  Storage, 

G.  F.  McMahon,  V.  R.  Bonner  and 
B.  S.  Eichert,  Feb  79,  32  pp. 

TP-61  Technical  Factors  in  Small  Hydropower 
Planning,  D.  W.  Davis,  Feb  79, 

35  pp. 

TP-62  Flood  Hydrograph  and  Peak  Flow  Frequency 

Analysis,  A.  D.  Feldman,  Mar  79  21  pp. 

TF-J3  HEC  Contribution  to  Reservoir  System 
Operation,  B.  S.  Bichert  and 
V.  R.  Bonner,  Aug  79,  28  pp. 

TP-64  Determining  Peak-Discharge  Frequencies  in 

an  Urbanizing  Watershed:  A  Case  Study, 

S.  F.  Daly  and  J.  C.  Peters,  Jul  79,  15  pp. 

TP-65  Feasibility  Analysis  in  Small  Hydropower 

Planning,  D.  W.  Davis  and  B.  W.  Smith, 

Aug  79,  20  pp. 

TP-66  Reservoir  Storage  Determination  by  Computer 
Simulation  of  Flood  Control  and 
Conservation  Systems,  B.  S.  Eichert, 

Oct  79,  10  pp. 

TP-67  Hydrologic  Land  Use  Classification  Using 
LANDSAT,  R.  J.  Cermak,  A.  D.  Feldman 
and  R.  P.  Webb,  Oct  79,  26  pp. 


ADA-106  2S8 


ADA-106  259 


ADA-109  757 


ADA-109  758 


ADA-109  759 


ADA-109  760 


ADA-109  761 


ADA-109  762 


ADA-109  763 


TP-68  Interactive  Nonstructural  Flood-Control 
Plannng,  D.  T.  Ford,  Jun  80,  12  pp. 


ADA-109  764 


TP-69  Critic*!  Water  Surface  by  Minimum  Specific 
Energy  Using  the  Parabolic  Method, 

B.  S.  Eichert,  1969,  IS  pp. 

TP-70  Corps  of  Engineers  Experience  with 
Automatic  Calibration  of  a 
Precipitation-Runoff  Model,  0.  T.  Ford, 

E.  C.  Morris,  and  A.  D.  Feldman, 

May  80,  12  pp. 

TP-71  Determination  of  Land  Us*  from  Satellite 

Imagery  for  Input  to  Hydrologic  Models, 

R.  P.  Webb,  R.  Cermak,  and  A.  D.  Feldman, 
Apr  80,  18  pp. 

TP-72  Application  of  the  Finite  Element  Method  to 
Vertically  Stratified  Hydrodynamic  Flow 
and  Water  Quality,  R.  C.  MacArthur  and 
W.  R.  Norton,  May  80,  12  pp. 

TP-73  Flood  Mitigation  Planning  Using  HEC-SAM, 

D.  W.  Davis,  Jun  80,  17  pp. 

TP-74  Hydrographs  by  Single  Linear  Reservoir 

Model,  J.  T.  Pederson,  J.  C.  Peters, 
and  O.  J.  Helweg,  May  80,  17  pp. 


TP-7S  HEC  Activities  in  Reservoir  Analysis, 

V.  R.  Bonner,  Jun  80,  10  pp. 

TP-76  Institutional  Support  of  Water  Resource 
Models,  J.  C.  Peters,  May  80,  23  pp. 

TP-77  Investigation  of  Soil  Conservation  Service 
Urban  Hydrology  Techniques, 

D.  6.  Altman,  W.  H.  Espey,  Jr.  and 
A.  D.  Feldman,  May  80,  14  pp. 

TP-78  Potential  for  Increasing  the  Output  of 
Existing  Hydroelectric  Plants, 

D.  W.  Davis  and  J.  J.  Buckley, 

Jun  81,  20  pp. 


ADA-951  599 


ADA-109  765 


ADA-109  766 


ADA-109  767 


ADA-109  756 


ADA-109  768 


ADA-109  769 


ADA-109  770 


ADA-109  771 


ADA- 109  772 


TECHNICAL  PAPERS  (TP) (Continued) 


HEC 

NUMBER 


TITLE 


TP-79 


TP-80 


TP-81 


TP-82 


TP-83 


TP-84 


TP-85 


TP-86 


TP-87 


TP-88 


TP-89 


Potential  Energy  and  Capacity  Gains  from 
Flood  Control  Storage  Reallocation 
at  Existing  U.  S.  Hydropower 
Reservoirs,  B.  S.  Eichert  and 

V.  R.  Bonner,  Jun  81,  18  pp. 

Use  of  Non-Sequential  Techniques  in  the 

Analysis  of  Power  Potential  at  Storage 
Projects,  G.  M.  Franc,  Jun  81,  18  pp. 

Data  Management  Systesm  for  Water  Resources 
Planning,  D.  W.  Davis,  Aug  81,  12  pp. 

The  New  HEC-1  Flood  Hydrograph  Package,  A.  D. 
Feldman,  P.  B.  Ely  and  D.  M.  Goldman , 
May  81,  28  pp. 

River  and  Reservoir  Systems  Water  Quality 
Modeling  Capability,  R.  G.  Willey, 

Apr  82,  IS  pp. 

Generalized  Real-Time  Flood  Control  System 
Model,  B.  s.  Eichert  and  A.  F.  Pabst, 
Apr  82,  18  pp. 

Operation  Policy  Analysis:  Sam  Rayburn 
Reservoir,  D.  T.  Ford,  R.  Garland 
and  C.  Sullivan,  Oct  81,  16  pp. 

Training  the  Practitioner:  The  Hydrologic 
Engineering  Center  Program, 

W.  K.  Johnson,  Oct  81,  20  pp. 

Documentation  Needs  for  Water  Resources 

Models,  W.  K.  Johnson,  Aug  82,  16  pp. 

Reservoir  System  Regulation  for  Water 
Quality  Control,  R.G.  Willey, 

Mar  83,  18  pp. 

A  Software  System  to  Aid  in  Making  Real-Time 
Water  Control  Decisions,  A.  F.  Pabst 
and  J.  C.  Peters,  Sep  83,  17  pp. 


ADA-109  787 


ADA-109  788 


ADA-114  650 


ADA-114  360 


AJA-114  192 


ADA-114  359 


ADA-123  526 


ADA-123  568 


ADA-123  558 


ADA-130  829 


ADA-138  616 


TECHNICAL  PAPERS  (TP) (Continued) 
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HEC  HBC  NTIS 

NUMBER  TITLE _  PRICE  NUMBER 


<2.00  Each 


TP-90  Calibration,  Verification  and  Application 
of  a  Two-Dimensional  Flow  Model, 

D.  M.  Gee,  Sep  83,  6  pp. 

TP-91  HEC  Software  Development  and  Support, 

B.  S.  Eichert,  Nov  83,  12  pp. 

TP-92  Hydrologic  Engineering  Center 
Planning  Models 
D.  T.  Ford  and  D.  W.  Davis, 

Dec  83,  17  pp. 

TP-93  Flood  Routing  Through  A  Flat,  Complex 

Floodplain  Using  A  One-Dimensional 
Unsteady  Flow  Computer  Program, 

J.  C.  Peters,  Dec  83,  8  pp. 


ADA-13S  668 


ADA-139  009 

ADA-139  010 


ADA-139  Oil 


TP-94  Dredged-Material  Disposal  Management 
Model,  D.  T.  Ford,  Jan  84,  18  pp. 


ADA-139  008 


TP-95  Infiltration  and  Soil  Moisture  Redistribution  ADA-141  626 

in  HEC-1 ,  A.  D.  Feldman,  Jan  84, 


TP-96  The  Hydrologic  Engineering  Center  Experience  ADA-141  860 

in  Nonstructural  Planning,  W.  K.  Johnson 
and  D.  U.  Davis,  Feb  84,  7  pp. 


TP-97  Prediction  of  the  Effects  of  a  Flood  Control  ADA-141  951 

Project  on  a  Meandering  Stream, 

D.  M.  Gee,  Mar  84,  12  pp. 


TP-98  Evolution  in  Computer  Programs  Causes  Evolution 
in  Training  Needs:  The  Hydrologic 
Engineering  Center  Experience,  V.  R.  Bonner, 
Jul  84,  20  pp. 


TP-99  Reservoir  System  Analysis  for  Water  Quality, 

J.  H.  Duke,  D.  J.  Smith  and  R.  G.  Willey, 
Aug  84,  27  pp. 


TECHNICAL  PAPERS  (TP) (Continued) 


HEC  HEC  NTIS 

NUMBER  TITLE _  PRICE  NUMBER 

>2.00  Each 

TP-100  Probable  Maximum  Flood  Estimation  -  Eastern 

United  States,  P.  B.  Ely  and  J.  C.  Peters, 

Jun  8A,  S  pp. 

TP-101  Use  of  Computer  Program  HEC-S  For  Water 
Supply  Analysis,  R.  J.  Hayes  and 
Bill  S.  Eichert,  Aug  84,  7  pp. 

TP-102  Role  of  Calibration  in  the  Application 
of  HEC-6 ,  D.  Michael  Gee,  Dec  84, 

19  pp. 

TP-103  Engineering  and  Economic  Considerations 
in  Formulating  Nonstructural  Plans, 

M.  W.  Burnham,  Jan  85,  16  pp. 

TP-104  Modeling  Water  Resources  Systems  for 
Water  Quality,  R.  G.  Willey, 

D.  J.  Smith  and  J.  H.  Duke, 

Feb  84',  10  pp. 

TP-105  Use  of  a  Two-Dimensional  Flow  Model  to 

Quantify  Aquatic  Habitat,  D.  Michael 
Gee  and  Daniel  B.  Wilcox,  Apr  85,  10  pp. 
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